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ABSTRACT 

Context. Geometrically thin, detached shells of gas have been found around a handful of carbon stars. The current knowledge on 
these shells is mostly based on CO radio line data. However, imaging in scattered stellar light adds important new information as well 
as allows studies of the dust shells. 

Aims. Previous observations of scattered stellar light in the circumstellar medium around the carbon star U Ant were taken through 
filters centred on the resonance lines of K and Na. These observations could not separate the scattering by dust and atoms. The aim of 
this paper is to remedy this situation. 

Methods. We have obtained polarization data on stellar light scattered in the circumstellar medium around U Ant through filters 
(-H ' which contain no strong lines, making it possible to differentiate between the two scattering agents. Kinematic, as well as spatial, 

information on the gas shells were obtained through high-resolution echelle spectrograph observations of the KI and NaD lines. 
Results. We confirm the existence of two detached shells around U Ant. The inner shell (at a radius of w 43" and a width of » 2") 
consists mainly of gas, while the outer shell (at a radius of x50" and a width of «7") appears to consist exclusively of dust. Both 
' shells appear to have an over-all spherical geometry. The gas shell mass is estimated to be 2 x 10~ 3 M , while the mass of the dust 

shell is estimated to be 5 x 10~ 5 M Q . The derived expansion velocity, from the KI and NaD lines, of the gas shell, 19.5 kms~', agrees 
with that obtained from CO radio line data. The inferred shell age is 2700 years. There is structure, e.g. in the form of arcs, inside the 
gas shell, but it is not clear whether these are due to additional shells. 

Conclusions. Our results support the hypothesis that the observed geometrically thin, detached shells around carbon stars are the 
results of brief periods of intense mass loss, probably associated with thermal pulses, and subsequent wind-wind interactions. The 
separation into a gas and a dust shell, with different widths, is most likely the effect of different dynamical evolutions of the two media 
after their ejection. 
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(N 
<N 

1. Introduction Olofsson et al. d 1988t found two stars with distinctly double- 

ts . peaked line shapes, suggesting from the star detached gas shells, 

O ■ In their final stages of evolution stars between » 0.8 - 8 M Q as- . m eflfect of isodic m ^ The s wafj extendedi 

cend the asymptotic giant branch (AGB). The evolution in this and a few additional objects were detected with sim0ar si 

phase is strongly affected by intense mass loss from the stel- of M M time . variable mass loss [summarize d in Olofsson et 

^ . lar surface with winds corresponding .to mass-crates in the d ^ i 996)] . Maps of the CO(/ = 1-0, 2-1, and 3-2) emis- 

H ; range 10~ 8 M yrj 1 to 10 4 M Q yr 1 (Blocker EES Although sion ^ revealed etrica „ ^ co line -emitting shells of 
W , the existence of the mass loss is well established much remains arQund ±Q carbon starfj R ^ {J ^ § § y644 s and 

in the understanding of the mechamsm(s) behind it. Present ob- ^ rr ., e , , i ln n<:t u- u i *■ a m, 

? . . , , TT Cyg (Olofsson et al. 1996). High resolution maps, made with 

servational uncertainties in the determined mass-loss ratesare the TRAM plateau de Bure inter fe r0 meter, of TT Cyg (Olofsson 

in the best case a factor of about three (Ramstedt et al. 2008), et al ffl and the carbon star v Cam (Lindqvist et al . 1999) 

based on smooth, spherical symmetric wind models. This is un- showed ^ ^ shdls ^ ^ {AR/R < Q l% and remarkabl 

fortunate, as even a moderate change in the mass-loss rate will herical with the detection of a CQ line .emitting shell around 

have a profound effect on the evolution of the star, its nude- the car bon star DR Ser, a total of seven cai'bon stai-s with geomet- 

osynthesis, and its return of chemically enrichedmatenal to the rfcj% ^ detached shells ^ knQwn (Sch5ier gt al 

mt ^ stella LEEd m T m (Forestall & Charbonnel ESI Schroder & Detached shells of dust have been observed around a number 

Sedlmayr 2001 ). In particular, the dependances of the mass loss of AGB afld t . AGB stars (Waters et al rjggj Izumiura et 

on stellar parameters such as mass, luminosity, radius, tempera- a , rjg^ rjggg s k et al M Warei et al These 

ture and pulsational properties are essentially unknown. dust ^ ^ ^ necessaril geometrically thin. R Hya is the 

In general, the mass loss is thought to increase(on aver- Qnl M AGfi star wkh a detected detached sheU> ^ Ms 

age) as the star evolves along the AGB (Habing[l996). In ad- case th h dust emission (Hashimoto et ai.rjg^. The 21-cm 

dition, temporal vanations on short as well as long time scales ffl Hne emission can fee used tQ stud mass . loss rate variations 

are present and they may be substantial. In a major survey of Qn a j ^ scale than ^ bed fe ±e CQ ^ emission 

circumstellar CO radio hne emission from nearby carbon stars, afj weU as effects of interaction with the ISM (e G6rard & Le 
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Bertre 2003). Also in this case shells are found, but they are of a 
different nature than the detached CO shells. 

A phenomenon which may affect the mass-loss properties 
during AGB evolution is the He-shell flash or thermal pulse. 
During a thermal pulse the star undergoes a structural change 
induced by explosive He-burning in the He-shell. The thermal 
pulse causes the surface luminosity, radius, and effective tem- 
perature of the star to change (Pols et al. 1200 It . The most sig- 
nificant features of the pulse are a rapid dip of the luminosity 
during the pulse, a very brief luminosity peak immediately after 
the pulse, followed by a slow dip and a subsequent long-lasting 
period when the star slowly recovers to its pre-pulse luminosity 
(Wagenhuber & Groenewegen 1 1 9981 ). The amplitude and dura- 
tion of the luminosity changes depend mainly on stellar mass 
and metallicity (Boothroyd & Sackmann 1 19881 1. The changes in 
the stellar structure mix regions with nuclear-processed material 
in the interior of the star with the surface through convection. 
This dredge-up leads to a chemical evolution of the stellar enve- 
lope, and may result in the formation of a carbon star, changing 
the C/O-ratio from < 1 (M-type AGB stars) to > 1 (carbon stars). 
Hence, an understanding of the thermal pulse cycle is essential 
for understanding AGB evolution. 

Olofsson et al. (1990) suggested that the geometrically thin, 
detached shells seen in CO line emission may be connected 
to an increase in mass-loss rate caused by the luminosity and 
temperature changes during a thermal pulse, a scenario devel- 
oped in more detail by Schroder et al. ( 11999b . However, Steffen 
et al. ( 119981 ) showed, using numerical models of the expand- 
ing matter, that an increase of mass-loss rate during the ther- 
mal pulse alone is not sufficient to form the observed detached 
gas shells. Instead, the interaction of a faster wind, associated 
with a brief increase in mass-loss rate, colliding with a previous 
slower wind can form geometrically thin shells of gas (Steffen 
et al. IT9981 Steffen & Schonberner |2"0"00j >. This is confirmed 
by detailed hydrodynamic models of the thermal pulse, the dy- 
namical stellar atmosphere, and the evolution of the expanding 
medium (Mattsson et al. 2007). A variation in mass-loss rate to- 
gether with a change in the expansion velocity is required to 
form the observed detached shells. Schoier et al. ( 120051 ) mod- 
elled the thermal dust and CO line emission for all known carbon 
stars with detached shells, and found an increase in shell mass 
and decrease in shell expansion velocity with increasing shell ra- 
dius, i.e., increasing shell age, indicating a two- wind interaction 
scenario where a faster wind sweeps up material from a previous 
slower wind. Interaction of the stellar wind with the surrounding 
matter (Libert et al. 2007) or bow shocks between the interstel- 
lar medium and the expanding circumstellar envelope (Wareing 
et al. 12006b may also lead to the formation of detached shells 
around some sources, although they are not geometrically thin. 

Although powerful in detecting the detached shells of gas 
around carbon stars, observations of CO radio line emission 
have their limitations. The emission depends on the excitation 
and chemistry of the CO molecules, and the translation from 
brightness distribution to density distribution is not straightfor- 
ward. Furthermore, the lack of interferometers in the southern 
hemisphere capable of observing the CO rotational lines limits 
the possibility to obtain detailed interferometer maps to the few 
northern sources. Infrared emission from dust provides a probe 
of the shells, however, these observations are usually of low spa- 
tial resolution. 

Gonzalez Delgado et al. ( T200T1 hereafter GD2001) used a 
novel technique to study the circumstellar medium around the 
carbon stars R Scl and U Ant. They imaged the circumstel- 
lar scattered stellar light through two narrow band (5nm) fil- 



ters (centred on the resonance lines of K and Na). By placing a 
choronograph over the star, the direct stellar light was avoided 
and the detached shells could be detected. This allowed for the 
first time to study the detailed structure of the shells taking ad- 
vantage of the high angular resolution obtained in optical data. 
The same technique was used again two years later, this time ob- 
serving the scattered stellar light from R Scl and U Ant in polari- 
sation mode (Gonzalez Delgado et al.|2Q031 hereafter GD2003). 
The major fraction of the polarised light is due to scattering by 
dust particles, and these observations gave high spatial resolu- 
tion images of the dust distribution in the circumstellar medium. 
In the case of U Ant they proposed a structure of four shells. 
The inner two shells were tentatively seen only in the images 
through the filter containing the NaD lines. For the outer two 
shells they found a difference in polarisation degree. Models of 
the polarised light indicated that the light from outermost shell 
is entirely due to dust scattering. However, the bulk of the total 
scattered intensity (a:70%) must come from another scattering 
agent, most likely due to atomic scattering (in the KI and NaD 
lines). 

The observations by GD2001 and GD2003 had the limitation 
that they could not conclusively separate the different contribu- 
tions by dust and line scattering. In this paper we present new 
observations of the circumstellar medium around U Ant taken 
with the EFOSC2 (ESO Faint Object Spectrograph and Camera) 
instrument on the ESO 3.6 m telescope through filters that con- 
tain no strong lines, hence making it possible to differentiate be- 
tween dust- and line-scattered light (new images in the filters 
centred on the resonance lines of K and Na were also obtained). 
In addition, we observed U Ant using the echelle spectrograph 
EMMI (ESO Multi-Mode Instrument) on the ESO NTT, and ob- 
tained a map of the circumstellar medium in CO(7 = 3-2) emis- 
sion with the APEX (Atacama Pathfinder Experiment) telescope. 
In Sect. [2] we present the basic properties of U Ant as well as 
describe the EFOSC2, EMMI, and APEX observations. We de- 
scribe the reduction process in Sect. [3j and present our results 
in Sect. [4] Finally, we discuss the properties of the circumstellar 
medium around U Ant and its origin in Sect. |5j where we also 
present our conclusions. 

2. Observations 

2.1. UAnt 

U Ant is an N-type carbon star with irregular variability. It is lo- 
cated at a distance of 260 pc (the Hipparcos distance; Knapp et 
al. |2003l >. It has a luminosity of 5800 L B , bas ed on SED fitting 
at wavelengths < lOyum (Schoier et al. I2005l l. The present-day 
mass-loss rate, based on CO radio line emission models, is low 
(« 2 x 10~ 8 M yr" 1 ; Schoier et al. [2003T >. Izumiura et al. ( fT997b 
estimate the ZAMS mass of U Ant to be 3 - 5 M Q , but this must 
be regarded as highly uncertain since it is based on the time dif- 
ference between two ejected dust shells assuming they are due 
to two consecutive thermal pulses. 

2.2. EFOSC2 imaging of polarised light 

The circumstellar environment of U Ant was observed in po- 
larised, scattered light during three nights in April 2002 us- 
ing the EFOSC2 focal reducer camera on the ESO 3.6 m tele- 
scope. The images were taken through an Ha filter (657.7 nm), 
a Strtimgren y filter (548.2 nm, hereafter Str-y), and narrow fil- 
ters (5 nm) centred on the resonance lines of Na (589.4 nm, 
hereafter F59) and K (769.9 nm, hereafter F77), with a pixel 
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Table 1. Observations of the circumstellar environment of U Ant. The table gives the instrument, the date of observation, the filter 
(or frequency) used, the filter width, the pixel size (for APEX the beam FWHM is given), the spectral resolution, and the total 
integration time. 



Instrument 


Date 


filter 


[nm] 


A^FWHM 

[nm] 


pixel 
["] 


Av 
[km s" 1 ] 


f 
[s] 


EFOSC2 


Apr. 2002 


Str-y 


548.2 


18.2 


0.32 




5000 






F59 


589.4 


5 


0.32 




4600 






Ho- 


657.7 


6.2 


0.32 




3900 






F77 


769.9 


5 


0.32 




90 


EMMI 


Mar. 2004 


F59 


589.4 


5 


0.33 


4.7 


9000 






F77 


769.9 


5 


0.33 


4.7 


9000 


APEX-2a 


Dec. 2006 




345 GHz 




18 


0.5 


15500 



a The total integration time for each polarisation angle is given for the EFOSC2 observations. For the APEX observations the total integration 
time for the map is given. 



scale of 0732/pixel. The average seeing during the three nights 
was » l'/3 (see Table[TJfor details on the observations and filters 
used). Since the stellar light outshines the circumstellar scattered 
light by a factor of » 10 4 , the use of a coronographic mask is 
necessary. A mask with a radius of ss 4" was chosen, reducing 
the direct stellar light enough to allow for long exposures with- 
out saturating the CCD, while still making it possible to detect 
scattered light close to the star. 

Polarimetric images were taken using the Ha, Str-y, and F59 
filters, using a rotating half-wave plate and a fixed polariser. 
Images were taken at polarisation angles of 0°, 45°, 90°, and 
135°. The total integration time was 3900s/angle, 4600s/angle, 
and 5000s/angle in the Ha, Str-y, and F59 filters, respectively. 
Standard stars were observed through all four filters (LTT3218 
in the Ha, F59, and F77 filters, and LTT6248 in the Str-y and 
F77 filters) in order to obtain the absolute flux calibration. In 
addition, a template star of similar magnitude and spectral type 
as U Ant (HD 137709) was observed in all filters, in order to 
characterise the stellar psf, in particular over the area of the cir- 
cumstellar envelope. An image in direct imaging mode without 
polarising filter was taken through the F77 filter. However, due to 
the limited available observing time, only one such image could 
be taken, with a total integration time of 90s. The resulting im- 
age is of rather poor quality, and we instead use the results in 
the F77 filter by GD2003. Finally, in order to determine the total 
stellar flux, U Ant was observed in direct imaging mode without 
the use of a coronograph in all filters. 



is dominated by a spherical, geometrically thin, expanding shell, 
this results in very characteristic elliptical shapes. The length of 
the observed ellipse gives a measure of the size of the circumstel- 
lar shell, while the width of the ellipse measures the expansion 
velocity of the shell (see Sect. 14. 5i . 

Spectra were taken at slit offsets of 15" east and west, and 
25" east of the star. The total integration time was 9000 s in 
each filter at each offset. The average seeing during the night 
was « 1". The pixel scale of EMMI is 0733/pixel. We used an 
078 wide slit, resulting in a spectral resolution of R ss 63500, or 
»4.7kms _1 . 

The standard star HR718 was unfortunately only observed 
during an observing run one year earlier when the same instru- 
mental set-up was used. This, in addition to the non-standard 
instrumental set-up of the observations, made a reliable flux- 
calibration impossible. 

2.4. APEX CO radio line emission 

An on-the-fly (OTF) map of the U Ant circumstellar envelope in 
the CO(7 = 3-2) line at 345 GHz was obtained with the APEX 
telescope equipped with the APEX-2a receiver during December 
2006. Due to the limited amount of observing time, we observed 
only the south-eastern quadrant in a 70" x 70" map. The total 
integration time of the map is 4.3 hours, resulting in an rms noise 
level of as 0. 1 K at a resolution of 0.5 km s . The beam width of 
APEX at 345 GHz is 18". 



2.3. EMMI echelle spectroscopy 

In March 2004 we performed spectroscopic observations of 
U Ant using the ESO NTT equipped with the echelle spectro- 
graph EMMI. Instead of using a cross disperser to separate the 
orders, we used the F59 and F77 filters to select the orders which 
contained the resonance lines. This allowed the use of a long slit 
that covered the entire shell. 

The slit was placed offset from the star, avoiding (most of) 
the stellar emission and making it possible to observe the cir- 
cumstellar scattered light instead. The expansion of the shell 
causes a splitting of the line. Light from the front part of the 
shell (as seen from Earth) is shifted to shorter wavelengths, while 
light from the rear part is shifted to longer wavelengths. The ve- 
locity difference between the two parts is largest where the line 
of sight goes through the middle of the shell, where the expan- 
sion along the line of sight is largest, while it is zero at the top 
and bottom edges of the shell. Since the circumstellar envelope 



3. Data reduction 

3.1. EFOSC2 images 

The EFOSC2 images were reduced using standard tasks for 
CCD reduction in IRAF (in particular the ccdproc task in the 
imred.ccdred package). The images were bias-subtracted, flat- 
fielded, and average-combined after aligning the individual im- 
ages using background stars. In order to subtract the stellar psf 
from the final images, it is important to know the exact position 
of the central star behind the occulting mask. The position was 
determined by finding the centre of rings of constant intensity 
close to, but outside of, the occulted region. The position was 
confirmed by measuring the distance to background stars in the 
image. 

The images were flux calibrated using the observed stan- 
dard stars. The surface brightness of the scattered light was 
calculated by determining the scale factor / sca ] e that converts 
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[counts s 1 pix 1 ] to [erg s 1 cm 2 ' 
culated as 



scale 



^std x AA 
^cnts X pixscale 2 ' 



]. The scale factor was cal- 



(1) 



where Fs t d is the tabulated flux density of the standard star in 
[erg s _1 cirT 2 A -1 ], AA the filter width in A, F cnts the number of 
counts per second in the image of the standard star (determined 
by a gaussian fit to the observations), and pixscale the pixel scale 
of the image in arcsecs. The flux calibration is estimated to be 
accurate within a factor of two. The subtraction of the template 
star, however, may lead to significant uncertainties in the deter- 
mined fluxes (see Sect. 13.21 ). The stellar fluxes were determined 
by directly comparing the number of counts per second in the 
standard star with the tabulated flux and using this factor to con- 
vert from counts per second to flux in the direct imaging frames 
of U Ant. 

3.2. Polarisation with EFOSC2 

Several different sources may contribute to the measured po- 
larised light: the central star, stellar light scattered in the inter- 
stellar medium, the sky, and the telescope and instruments, as 
well as polarised light from the sky background. In the follow- 
ing analysis, however, we assume that the stellar light is unpo- 
larised. This is confirmed by BVRI polarimetry studies of vari- 
able stars, setting a strong upper limit on the polarisation of the 
stars (< 1%, Raveendran 1 1 99 1 \ . Further, the relative proximity 
of U Ant suggests a negligible contribution to the observed po- 
larisation from the interstellar medium. The observations were 
made just after new moon, limiting the amount of polarised light 
due to the background. Finally, although time limitations pre- 
vented us from taking images of polarimetric standards in order 
to characterise the effects of the telescope and the instrument, 
there are no indications that these might introduce any signif- 
icant polarisation in the observations. This was confirmed by 
GD2003, who observed template stars using the same instru- 
mental set-up. They found no polarisation in the regions where 
the scattered stellar light was detected. 

Therefore, assuming that all the polarised light is due to scat- 
tering in the circumstellar environment of the star, the Stokes 
parameters for the shell emission are given by 



hh ~ 

Qsh 

Psh '■ 



Is, - I, 



■ Qm 

■ u m 



(2) 



where /„,, /„, and 1^ are the total measured intensity, the total 
intensity due to stellar light scattered in the interstellar medium, 
Earth's atmosphere and the instrument, and the total intensity in 
the sky background, respectively. Q m and U m are the measured 
Stokes' parameters, and P s i, is the total polarised flux in the shell. 
The intrinsic polarisation degree of the shell is given by 



Psh - 



Li i 



(3) 



The Stokes' parameters can now be determined in a straight- 
forward way from the observed images 



I m = F(0°) + F(90°) 
Q m = F(0°) - F(90°) 
U m = F(45°)-F(135°), 



(4) 




10 20 30 40 50 60 70 80 
Distance to the star [arcsec] 

Fig. 1. Azimuthally averaged radial profiles (AARPs) in the F59 
filter of U Ant (solid line), the template star (dashed line), and a 
moffat function fitted to the data (dotted line). 



where F indicates the frames taken at the different polarisation 
angles. 

The most difficult task is the determination of the total in- 
tensity in the circumstellar shell /.,/,. The stellar psf, due to light 
scattered in the Earth's atmosphere and the telescope, was first 
removed by scaling and subtracting an image of a template star 
with the same spectral characteristics as U Ant. This requires 
an exact alignment of the template and source stars, as well as 
a very well-defined and 'clean' template psf. Unfortunately, the 
images of the template star contained a large number of back- 
and foreground stars, making it impossible to cleanly subtract 
the stellar psf. We therefore fitted Moffat profiles (Moffat [T9691 
to the azimuthally averaged radial profiles (AARPs) of the 
brightness distributions in the U Ant template star images (the 
latter were averaged over the same position angles as in the 
U Ant images). The Moffat profile is described by 



/ = i x [i + {-fr b , 

ro 



(5) 



where Iq is the maximum intensity, r the distance from the cen- 
tre, and ro the width of the profile. The exponent b determines 
how quickly the profile decreases. The psf is assumed to be 
circularly symmetric. The best-fit Moffat profile is determined 
manually taking both the U Ant and the template star profiles 
into consideration. The exponent, width, and maximum intensity 
of the Moffat profile were adjusted until it matched the U Ant ob- 
servations close to the edge of the occulting mask (< 10", where 
the stellar flux dominates over the circumstellar scattered light), 
and the tail of the psf well outside of the circumstellar region 
{<: 70"). Figure Q] shows the AARPs of the total observed inten- 
sity in U Ant and the template star, and the fitted Moffat profile 
in the F59 filter. The Moffat profile fit to the template star works 
reasonably well, and the excess emission due to the circumstel- 
lar envelope in the U Ant data can clearly be seen. However, the 
fit to the U Ant data is not as good at large distances from the 
star, and the shape of the AARP of the total shell intensity is very 
sensitive to the fit of the moffat function. 

The uncertainty of the template subtraction and other instru- 
mental effects significantly reduce the quality of the information 
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in the central parts of the images. As reliable measurements are 
not possible, this region is blanked out in the images by a mask 
with a 20"radius. The dominant diffraction spikes due to the spi- 
ders were blanked out as well. 

3.3. EMMI echelle spectra 

The spectroscopic images from EMMI were reduced using stan- 
dard tasks in IRAF for flat fielding and bias subtraction. The 
spectrum obtained on the CCD is usually significantly curved 
when a long slit is used. This curvature is not necessarily con- 
stant over the entire CCD, and the spectrum needs to be carefully 
straightened. This was done using different tasks for long-slit 
spectroscopy in the IRAF noao.twospec.longslit package. 

Due to the use of order-sorting filters instead of a cross dis- 
perse^ the determination of the dispersion in the EMMI data 
proved to be problematic. In fact, three lines were visible in the 
F77 data, although at most two (the two resonance lines) were 
expected, indicating an order overlap. This was confirmed by 
examining the profiles of the flat-field images along the disper- 
sion axis, showing that the F77 filter has its centre in the region 
between two orders. Hence, the lines observed in the F77 fil- 
ter show two orders of the 769.9 nm line, and one order of the 
766.5 nm line. Although the dispersion usually does not vary 
much between neighbouring orders in echelle spectra, the widths 
of the ellipses (due to the expansion of the shell) in the two or- 
ders differs significantly (by a 25%). The widths of the 769.9 
nm and 766.5 nm lines in the same order differ as well, indi- 
cating that the dispersion varies also within one order. The two 
orders lie close to the edge of the CCD, hence the dispersion 
may be affected by the optics of the instrumental setup. A care- 
ful determination of the dispersion at the position of the lines is 
therefore extremely important for a correct determination of the 
kinematical properties. 

The identification of spectral lines in the calibration spectra 
is nearly impossible when different orders overlap. However, al- 
though the images were taken well offset from the central star, 
stellar light scattered in the atmosphere and the instrument is still 
projected onto the slit. This, in principle unfortunate fact, made 
it possible to compare the stellar spectrum with an archived and 
wavelength-calibrated spectrum of a star of similar spectral type 
(the carbon star HD20234), resulting in a good dispersion cali- 
bration. 

The centre of the F59 filter lies close to the centre of one or- 
der, hence there is little or no order overlap in this filter. In this 
case, however, the calibration spectrum is dominated by one very 
strong line, making it impossible to identify other lines in the 
spectrum for a wavelength calibration, fin the case of ordinary 
echelle spectra, this does not create a problem, as the dispersion 
solution is extrapolated from the calibration in neighbouring or- 
ders.] The atmospheric Na resonance lines are clearly visible in 
the data though, and the comparison with a spectrum of a star of 
similar spectral type leads to a good dispersion solution also in 
this filter. 

3.4. APEX CO radio line data 

The on-the-fly data was reduced using CLAsfl The entire data 
set was inspected manually, and about 80 spectra out of 735 
were dropped because of various anomalies. The baselines in 

1 Continuum and Line Analysis Single-dish Software from 
the Observatoire de Grenoble and Institut de Radio-Astronomie 
Millimetrique 




RA offset ["] 

Fig. 2. A CO(J = 3-2) velocity-integrated intensity map in a 
+ 5kms _1 interval centered on the systemic velocity. U Ant is 
located in the (0,0)-position and the map covers the second quad- 
rant of the detached shell. Emission from the detached shell as 
well as the present-day mass-loss wind is apparent. The contours 
are based on a grid-spacing of 8", while the colour image has an 
interpolated grid spacing of 4". The contour levels are 0.3, 0.6, 
0.9, 1.2, and 1.4 K km s -1 . 

the remaining data were relatively stable, and only second order 
polynomial fits were subtracted. Once we had a consistent data 
set, the remaining spectra were averaged per position, using a 
weight determined by the noise level, to a grid-spacing of 8". 
Figure [2] shows the map of the velocity-integrated intensity, in a 
+ 5 km s _1 interval centered on the systemic velocity (this range 
is chosen since it gives the best view of the shell position, con- 
sidering also the signal-to-noise ratio), in the second quadrant 
of the detached shell. The shell can clearly be seen and its lo- 
cation is consistent with the results from the CO{J =1-0 and 
2-1) maps reported by Olofsson et al. (1 19961 1. Emission from 
the present-day mass loss is also visible. The decrease of emis- 
sion to the north and west is an artefact due to the limited map 
size in combination with the beam size. 

4. Results 

4.1. Review of previous results 

CO radio line observations of U Ant show a, not so common, 
highly double-peaked profile. This was interpreted as emission 
from a geometrically thin, detached gas shell around U Ant, 
an interpretation confirmed by mapping the CO line emission 
(Olofsson et al. ll996t . Models show that the CO data is consis- 
tent with a detached shell of gas with a radius of « 43" (Schoier 
et al. [20051 . GD2001 and GD2003 observed the circumstellar 
environment of U Ant in both direct imaging and polarisation 
mode through two filters covering the Na and K resonance lines 
(the F59 and F77 filters, respectively). They introduced four 
shells, at « 25", 37", 43", and 46" (shells 1 to 4, respectively). 
Shells 1 and 2 were only tentative. Shell 4 was introduced ten- 
tatively in GD2001, and it was confirmed in GD2003. The po- 
larisation measurements of GD2003 indicate a separation of the 
dust and the gas, shell 4 consisting of mainly dust while shell 3 is 
dominated by gas. The CO radio emission lines most likely come 
from shell 3. Figure [3] shows the positions of the different shells. 
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Fig. 3. A schematic drawing showing the positions of shells 
around U Ant as derived in GD2001 and GD2003 (shells 1 to 
4 in order of increasing size; shells 1 and 2 were tentatively in- 
troduced). The vertical lines indicate the positions of the slit in 
the EMMI observations presented here. 



The vertical lines show the positions of the slit in the EMMI 
data presented here. Izumiura et al. ( 119971 ) reported the detec- 
tion of two dust shells in IRAS images at 60 fi and 100 fim. They 
processed the images to obtain higher-resolution IRAS (HIRAS) 
images and derived shell radii of ss 46" and ss 3', respectively. 
The shell at 3' is outside our field of view, and will not be con- 
sidered further. 

In light of these previous results, we present the results of 
our new data in the following sections, describing the images 
in polarised, scattered light from EFOSC2 (Sect. 14. 2\ and the 
derivation of shell fluxes, radii and widths (Sect. 14.4b . The ex- 
pansion velocity of shell 3 is determined from the EMMI data 
(Sect. [43]). The different contributions to the scattered light from 
the dust and gas are estimated (Sect. 14.6b and the new CO data 
from APEX is modelled (Sect. 14.71 ). Finally, the dust masses in 
the individual shells are determined (Sect. |4~8l l. 

4.2. Images and radial profiles of the scattered light 

Figure H] shows the total flux, I S h, and polarised flux, P s h, images 
in the Str-y, F59, and Ha filters. AARPs of the I sh , P, h , and the 
polarisation degree, /?,/,, are shown in Fig. [5] 

The morphology in the Ha and Str-y I s u images differs sig- 
nificantly from that in the F59 filter. The former both show a cir- 
cular disk of ss 50" radius. The I s h AARPs in both filters show 
that the intensity remains relatively constant out to a radius of 
~ 50" followed by a tail of declining intensity. GD2003 showed 
that such a brightness distribution is consistent with scattering 
of stellar light in a detached shell of dust, assuming that the dust 
has a preferential scattering efficiency in the forward direction. 
The irregularities in the AARPs are likely due to a few bright 
arcs present in the images closer to the star. However, a small 
peak at = 43" indicates the presence of an additional shell, and 
an extended arc at this distance from the star is clearly seen in 
both filters in the south-east quadrant. 



The / s h AARP of the F59 image shows a much more pro- 
nounced peak at ss 43", while the tail at ss 48"is relatively weak. 
Since the F59 filter may contain a contribution from stellar light 
scattered in the NaD resonance lines, while the Ha and Str-y 
filters are dominated by dust scattered light, this confirms the re- 
sults by GD2003, where shell 3 is dominated by line scattering, 
while shell 4 is dominated by dust scattering. 

Only light that is scattered at an angle of ss90° will be 
strongly polarised. The polarisation due to the circumstellar 
medium will hence be strongest in the plane of the sky that goes 
through the star. Detached, spherical shells appear as ring-like 
structures in such images, directly revealing the spatial struc- 
ture of the shells. The images in the polarised flux P s j, in all 
filters clearly show a nearly perfect circular geometry. The cor- 
responding AARPs show the presence of two shells at ss43" 
and ss 48", corresponding to the positions of shells 3 and 4 in 
GD2003, but GD2003 did not detect shell 3 in the polarised light 
images. These AARPs make it possible to determine the loca- 
tions and the widths of shells 3 and 4. Fits to the AARPs of the 
P s h images in all filters show that the total polarised flux is dom- 
inated by shell 4 (see Sect. 14.41 ). Shell 4 also clearly dominates 
the AARPs of the p s f, images, the degree of polarisation reach- 
ing a maximum of 25-30% in all three filters at a radius of «50", 
confirming that shell 4 is dominated by stellar light scattered by 
dust. Since p s i, is derived P s h/hh> the outer parts (where I s i, and 
P s h ss 0) become less reliable. Hence, the picture of a shell of 
gas at the position of shell 3 containing a small amount of dust, 
and a detached shell dominated by dust with almost no gas at the 
position of shell 4 is strengthened. 

Structure inside shell 3 is apparent in the F59 image (but not 
in the others, including the F77 image of GD2003). The posi- 
tions of the peaks in the AARP are consistent with shells 1 and 
2 in GD2001 at 25" and 37". It is not clear, however, whether 
the observed structure is due to clumpy structures in shells 3 
and/or 4, or due to additional detached shells closer to the star, 
see Sect. 15. 1.11 

4.3. Scattering in detached gas and dust shells 

In order to understand the origin of detached shells, it is impor- 
tant to accurately determine their physical parameters, such as 
e.g. shell radius and width, and its density distribution and possi- 
ble dumpiness. Any determination of physical properties of the 
shells based on the data presented here requires the proper treat- 
ment of the scattering agent properties. For an individual dust 
grain these are affected by, e.g., the grain composition, size, and 
shape. The observed brightness distributions (due to scattering 
by a large number of grains) are affected by the grain size distri- 
bution, the density distribution within the shell, and the size and 
width of the shell. 

The scattering by the grains can be calculated using Mie the- 
ory, given the optical constants of the grains. If the grain size 
is small compared to the scattering wavelength (a/A <sc 0.16), 
the condition for Rayleigh scattering is fulfilled and the scatter- 
ing can be assumed to be isotropic. However, typical grain sizes 
are « 0.1 /jm, i.e. comparable to the wavelengths in our data, and 
asymmetric scattering becomes important. In order to determine 
the amount of scattering in a particular direction (i.e. the proba- 
bility for scattering in a direction with respect to the forward 
direction), we use an analytical expression (based on light scat- 
tered in the interstellar medium; Henyey & Greenstein |1941| ) 

m = rTT 2~/ M2 - (6) 

(1 + g z - 2gcos6y iz 
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Fig. 4. The results of the EFOSC2 observations of U Ant. Top row, left to right: the total intensity (/,/,) in the Str-y, F59, and Ha 
filters, respectively. Bottom row, left to right: the polarised intensity (P s i,) in the Str-y, F59, and Ha filters, respectively. 



where g is the scattering asymmetry parameter defined as 
£l(6) cosO sine dO 

£ i(e) sme de 



(7) 



where 1(6) is the scattering phase function. For g — the scatter- 
ing is isotropic, while for g — 1 all light is scattered in the forward 
direction and for g — -1 all light is scattered backwards. Draine 
( 120031 1 re-examined the expression by Henyey & Greenstein and 
concluded that it is valid at wavelengths between 0.4 /an and 
1.0 /mi, i.e., the wavelength range of interest here. The amount 
of polarisation as a function of the scattering angle is given by 



1 - cos 2 6 
1 + cos 2 £ 



(8) 



The left panel of Fig. [6] shows the scattering efficiency Q scdt 
of dust grains and the scattering parameter g vs. wavelength for 
amorphous carbon grains (Suh 2000) assuming spherical grains 
with radii of 0.1 yum. Equation [6] can then be used to derive ra- 
dial brightness distributions for different wavelengths assuming 
single scattering. The middle panel of Fig. [6] shows the expected 
brightness distributions for optically thin scattering by dust in a 
detached shell located at R — 48" from the star. They are given 
for complete isotropic scattering (g = 0) and for the calculated g 
for the respective filters. The density distribution of the detached 
shell along the radial direction is assumed to follow a Gaussian 
distribution with a FWHM of AR = 7" peaking at R. The grains 



have a constant size of 0.1 pm. The right panel of Fig. [6] shows 
the corresponding profiles for the polarised light. 



For all cases the limb brightening due to the narrow shell is 
apparent at a 48 "(hereafter referred to as the 'peak'). However, 
as forward scattering becomes more effective for shorter wave- 
lengths, the brightness profile at line of sights closer to the star 
increases compared to the peak intensity. Note that the amount 
of forward scattering does not have an effect on the position of 
the peak or the sharp decrease at larger radii. However, the den- 
sity distribution within the shell affects the shape of the peak. A 
sharp inner edge moves the position of the peak somewhat closer 
to the star, while a sharp outer edge leads to a sharper decline at 
large radii. The decline in brightness at large radii follows to 
good approximation a Gaussian distribution and it gives reason- 
able estimates of the radius and width of the shell (Mauron & 
Huggins 2000). However, fitting a Gaussian distribution to the 
peak and the tail results in a slightly underestimated radius R 
(by *8%) and an overestimated FWHM of the shell (by «20%). 
Hence, the detailed grain size distribution and the density distri- 
bution across the shell are the main uncertainties in determining 
the shell radius and width from the observed brightness distribu- 
tions. 



In the case of line scattering we use Equ [6] with g = for 
calculating the brightness distribution. 
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Fig. 5. AARPs from the EFOSC2 observations of U Ant. Top row, left to right: the total intensity (7 S /,) in the Str-y, F59, and Ha 
filters, respectively. Middle row, left to right: the polarised intensity (P s /,) in the Str-y, F59, and Ha filters, respectively. Bottom row, 
left to right: The polarisation degree {p s h) in the Str-y, F59, and Ha filters, respectively. The lines are: the observed data (solid), fits 
to the individual shells (dotted), and the sum of the individual fits (dashed). See text for details. 



4.4. Decomposition of the AARPs 

In order to derive quantitative results for the individual shells it is 
necessary to decompose the AARPs obtained from the observed 
images. The images in polarised light best show the spatial struc- 
ture of the dust, and we use Equ. [6] and [8] to derive theoretical 
AARPs for the polarised scattered light. The scattering asym- 
metry parameter g is determined for each filter using Mie theory 
and assuming spherical, 0.1 /mi sized carbon grains (Suh 120001 . 
A Gaussian density distribution of the dust as described above is 
assumed. The model distributions are fit to the peaks in the 



AARPs, hence determining the radius and widths of shells 3 and 
4. Equ.|6]is then used to calculate model profiles of the total in- 
tensity 7 s h. In the F59 filter shell 3 is assumed to be dominated by 
scattering in lines, and hence isotropic scattering is assumed in 
this filter (g = 0). The dust contribution to 7 s h from shell 3 in this 
filter is neglected. The dotted lines in Fig. [5] show the individ- 
ual brightness distributions. Tables [2] and [3] show the results of 
fitting the brightness distributions to the polarised and total flux 
AARPs, respectively. For the latter, the shell radii and sizes de- 
termined from the polarisation data are used. Hence, we confirm 
the detection of detached shells at 43"( shell 3) and 50"(shell 4) 
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from the star. The width of shell 3 is = 2'.'2, while shell 4 is a fac- 
tor of 3 broader. The results for the F77 filter in Table [3] are taken 
from GD2003. The fluxes for the individual shells in Table[3]are 
derived by summing the model brightness distributions over all 
angles. The total fluxes for all shells are calculated by summing 
the observed profiles over all angles (labelled 'total' in Table[3]). 

In all measurements probing the contribution to the flux from 
dust scattered light (i.e. the polarisation measurements and the 
7 s h images in the Str-y and Ha filters), shell 4 clearly dominates 
over shell 3. The contribution to the total flux from shell 3 is only 
comparable to that from shell 4 in the 7 s h images in the F59 filter. 
This is likely due to the added contribution of line scattered light 
in shell 3 in this filter. 

The theoretical brightness distributions do not perfectly fit 
the observed AARPs. Artefacts in the images (e.g., diffraction 
spikes), the uncertainty due to the psf subnaction, and a pos- 
sible clumpy snucture in the detached shells are likely to af- 
fect the AARPs, in particular in the inner parts of the images. 
The assumption of a constant grain-size in the model AARPs 
and a Gaussian density distribution across the shells also affects 
the model brightness distributions. The gradual decline at larger 
radii could also be a kinematic effect, e.g., due to grains of differ- 
ent size travelling at different velocities. The fit to the I s h AARP 
in the F59 filter is particularly poor. This is partly due to the 
presence of structures inside shell 3, but there is most likely also 
an effect of the NaD line scattering being, at least partially, op- 
tically thick (see below). Also, the FWHM of shell 3 in the 7 sh 
AARP seems to be wider than the one determined from the 
AARP. This is likely due to the small amount of dust present in 
shell 3. 



4.5. Kinematic and spatial information from the EMMI data 

Examples of the long-slit NaD and KI line data taken with 
EMMI are shown in Figs Q and [8] respectively. The ellipse-like 
shape of the lines due to the expansion of a geometrically, es- 
sentially spherical thin circumstellar shell can clearly be seen. 
As mentioned above, the NaD lines show a more complicated 
morphology, revealing additional structure closer to the star. 
However, we focus the analysis on the dominating shell. Ellipses 
were fit to the data (obtained by making cuts along the disper- 



sion axis and fitting Gaussian profiles to the line intensity distri- 
bution) using the conical representation of an ellipse 



A-x 2 + B-xy + C-y z +D-x + E-y = l, 



(9) 



where x and y are the points of the ellipse along the spatial and 
dispersion axes, respectively. The parameters A, B, C, D, and 
E were determined by fitting the the ellipse to the centres of the 
Gaussian profiles using a least-squares-fit method. The length r e \ 
and width v e i of the fitted ellipse are then given by (in arcseconds 
and kmr 1 , respectively) 



.ID 2 E 2 
r e [ = pixscale x I — + ^-p- + 



4ACI 



0.5 



c-d 1 

Vel = — X I - + 



,1 



c 



D 2 E 2 
4AC + 4C 2 



0.5 



(10) 



(ID 



where c is the speed of light, d the dispersion in A pix~ l , and 
A the wavelength of the line. The ellipse parameters are cor- 
rected for the projection and hence converted to the size R and 
expansion velocity v exp of the shell. In the case of NaD, the dou- 
blet lines were combined in order to increase the signal-to-noise 
ratio. The F77 filter is wide enough to also detect the line at 
766.5 nm, but it is (significantly) weaker than the 769.9 nm line. 
In addition to the radius of the shell, a cut through the data along 
the spatial axis at the systemic velocity gives an estimate of the 
width of the shell AR. The resulting parameters from the dif- 
ferent ellipse fits and the resulting averages for the shell size, 
shell width (the shell width determined in the EMMI data is the 
FWHM of the line along the spatial axis), and expansion velocity 
are given in Table [4] 

The EMMI data is dominated by scattering in a shell with 
a radius of »40" and a width of a 3", that expands at a veloc- 
ity of 19.5kms~'. Although the fit of an ellipse to the data is 
generally very good, effects due to the physical width of the slit, 
and the resolved spatial structure of the shell result in deviations 
from a perfect ellipse at the end points. This causes the mea- 
sured size of the shell to be systematically smaller (by a 2") than 
suggested by the peaks measured directly in the data. Taking 
this into consideration, the results are fully consistent with the 
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Table 2. The results of decomposing the P s i, AARP data into shell brightness distributions. R and AR are the radius and FWHM of 
the Gaussian density distribution assumed for the shell, / pea k and F tot are peak intensity and total flux of the brightness distributions 
fitted to the data. Pol. deg. is the peak polarisation degree at the position of the shells. 



Filter 


shell 


R 


AR 


'peak 


^tot 


pol. deg. 






["] 


["] 


[10- 16 ergs-» cm- 2 "~ 2 ] 


[10 -13 ergs" 1 cm" 2 ] 


[%] 


Stry 


3 


43.6 


1.9 


1.6 


2.4 


23 




4 


49.6 


5.4 


3.7 


12.0 


26 


F59 


3 


43.4 


2.4 


0.3 


0.5 


8 




4 


49.6 


7.1 


0.9 


3.4 


24 


Ha 


3 


43.4 


2.4 


0.5 


0.8 


25 




4 


49.8 


7.1 


2.0 


7.1 


33 



Table 3. The results of decomposing the I s h AARP data into shell brightness distributions, /peak, F to u and 5 av are the peak intensity, 
total flux, and average flux density of the shell brightness distribution fits (see text for details). The values for the total of all shells 
are obtained by integrating the observed brightness distributions. The numbers within brackets in the F59 and F77 filters in the flux 
column give the percentage of flux that may come from dust, based on a linear inter- and extrapolation between the Str-y and Ha 
filters. Fs C /F+ gives the ratio between the circumstellar and stellar fluxes. 



Filter 


shell 


/peak 


-ftot 










[10- 16 ergs-' cm- 2 "- 2 ] 


[10- 13 ergs-' cm" 2 ] 


[10- 14 ergs-» cm" 2 A" 1 ] 


[10- 4 ] 


Stry 


3 


5.7 


16.1 


0.9 


4.1 




4 


14.2 


80.1 


4.4 


20.1 




total 




144.4 


7.9 


36.3 


F59 


3 


8.7 


20.5 (19%) 


3.9 


20.5 




4 


4.0 


24.2 (83%) 


4.6 


24.2 




total 




84.7 (43%) 


16.3 


84.7 


Ha 


3 


2.4 


6.6 


0.5 


1.6 




4 


6.1 


34.9 


2.8 


8.3 




total 




68.2 


5.5 


16.3 


F77 


3 




14.0 (3%) 


2.8 


3.3 




4 




3.3(100%) 


0.7 


0.8 




total 




30.0 (50%) 


6.0 


7.1 



a The fluxes of the star in the filters are: 4.0 x lO -9 erg s~' cnr 2 (Str y), 1.1 x lO -9 erg s~' cnr 2 (F59), 4.2 x 10 -9 erg s -1 cnr 2 (Ha), and 4.2 x 
10" 9 ergs" 1 cm" 2 (F77). 



vations also contain stellar light scattered by dust, and it is dif- 
ficult to disentangle the contributions from the different scatter- 
ing agents. Observations of polarised light mainly measure light 
scattered by dust (resonance line scattering only accounting for 
up to « 5% of the polarised light; Loskutov & Ivanov 2007)>, and 
observations with filters that contain no lines make it possible to 
determine the amount of dust-scattered light, and, together with 
observations in the "line"-filters, the amount of line-scattered 
light. 



Table [3] gives the flux densities in the different filters. 
Through linear inter- and extrapolatation of the data in the Ha 
and Str-y filters (assumed to contain contributions only from 
dust), it is possible to estimate the fractions of flux in the F59 
and F77 filters that are due to dust. The result is that «43% and 
~ 50% of the total fluxes in the F59 and F77 filters, respectively, 
can be attributed to scattering in dust. The total polarised flux is 
« 7 times higher in shell 4 than in shell 3 in all filters (Table 
showing that the dominant contribution to the dust scattered light 
must come from shell 4. This is in agreement with the results in 
the polarisation degree, reaching its maximum (« 30%) at the 
position of shell 4. Hence, we can conclude that shell 3 is dom- 
inated by gas (and line scattering clearly dominates in the F59 
(« 81% of the total flux) and F77 (* 97% of the total flux) filter 
images, while dust dominates in shell 4. 



EFOSC2 and CO radio line data on shell 3. The expansion ve- 
locity determined from the resonance line data is in excellent 
agreement with that obtained from models of the CO emission 
lines, 19.0 km s 1 (Schoier et al. |2003] >. 

The brightness of the red-shifted side of the shell (i.e., the 
rear part as seen from Earth) is clearly brighter than the front in 
the NaD data. The same effect can be seen to a lesser extent in 
the KI data. This most likely is an optical depth effect. A higher 
optical depth in the NaD lines leads to increased reflection of the 
light on the red-shifted side of the shell, while more of the light 
is scattered away from the line of sight on the blue-shifted side. 

Finally, the EMMI data taken through the F59 filter also 
shows structures inside shell 3 in the NaD line, but ellipses are 
not apparent (see Fig.[7]l. A cut along the spatial axis at the sys- 
temic velocity reveals peaks also at 23" and 34" (corrected for 
projection effects), at the positions of the tentative shells 1 and 
2. This structure is not seen in the F77 filter data, i.e., in the KI 
line. 

4.6. Line vs. dust scattering 

Previous observations of scattered light in the detached shells 
around U Ant (GD2001 and GD2003) were made in filters con- 
taining strong resonance lines (except for a tentative detection 
of shell 4 in a Stromgren b filter at 469.0 nm). Since the widths 
of the lines are much smaller than the filter width, such obser- 



M. Maercker et al.: The detached shells around U Ant 



11 



) T 




' .. "TT— ■.. -}.\ 


40 20 


-20 -40 


) 




V- . . .1 


40 20 


-20 -40 


)- 




, ; .. . ^ ( . 



40 20 -20 

offset [arcsec] 



-40 



Fig. 7. EMMI long-slit spectra towards U Ant showing the NaD 
resonance lines. The doublet lines are combined in order to in- 
crease the signal-to-noise ratio. The slit is offset from the star 
by 15"east (top), 15"west (middle), and 25"east (bottom). The 
velocity scale is given with respect to the systemic velocity. 



Table 4. The results of the EMMI data. The NaD doublet lines 
were added and averaged before fitting an ellipse (see text for 
details). <x s and cr v , exp are the uncertainties in the radius R and 
the expansion velocity v exp of the shell, repectively. AR is the 
shell width (measured as the FWHM of the line along the spatial 
axis). The results are corrected for the projection. 
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Fig. 8. Same as Fig. |7]but for the KI resonance line at 769.9 nm. 



4.7. The CO shell 

Compared to the CO(/ = 1-0 and 2-1) data in Olofsson et 
al. (11996b , the APEX data presented here provides a higher an- 
gular resolution (18" beamwidth). Figure [9] shows the inten- 
sity (averaged over the central lOkms -1 ) AARP of the new 
CO(J = 3-2) data. The emissions from the shell and the present- 
day mass-loss wind are clearly separated. We here model the 
new CO(J = 3-2) data to show that the CO shell coincides with 
shell 3. The models are based on the best-fit model presented by 
Schoier et al. d20G51 ), slightly adjusted to fit the new APEX data. 
The code used in the CO line modelling is described in Schoier 
& Olofsson d200~n >. 



Following Schoier et al. (120051 1 we adopt a shell width of 
1.6xl0 16 cm (this corresponds to 2'.'6 at the distance of U Ant, 
i.e., close to the shell width as estimated from the line scattering 
data, which is much smaller than the angular resolution of the 
APEX CO line data). The kinetic temperature in the shell is set 
to 350 K. The density in the shell is assumed to be constant and 
defined by the gas shell mass (see below). The expansion veloc- 
ity of the shell is estimated to be 19.0 km s~' by fitting the CO 
line profiles. For the present-day mass loss we use a mass-loss 
rate of 1 .2x 10~ 8 M Q yr~' . This is about a factor of two lower than 
that obtained by Schoier et al. ( 120051 ). The latter is a result of an 
average fit to three CO lines, and the SEST CO(7 = 3-2) line is 
almost a factor of two stronger than the corresponding APEX 
line. However, it is clear that the APEX line is more reliably cal- 
ibrated than the SEST line (SEST had a main beam efficiency as 
low as 25% at 345 GHz), and we therefore adjust the mass-loss 
rate to fit the APEX line. In fact, as can be seen both in Figs [2] 
and|9] the present-day mass-loss emission is well separated from 
the detached shell emission, and only contributes to the AARP in 
the form of a Gaussian profile with the same width as the APEX 
beam. Therefore, the main goal of the APEX data to verify the 
position of the detached CO shell is not dependent on the as- 
sumption of the present-day mass-loss rate. The stellar radiation 
field is a central blackbody with a luminosity of 5800 L and an 
effective temperature of 2800 K. Finally, we adjust the position 
of the shell to get the best fit to the CO(7 = 3-2) intensity AARP. 
The resulting shell radius is 41". 

The result of the modelling is presented in the form of a 
model AARP in Fig. [9] The observed AARP is consistent with a 
model of a detached shell with a radius of 41", i.e., the CO shell 
coincides, well within the uncertainties, with shell 3, the shell 
responsible for the atomic line scattering. The H2 mass of shell 
3 is estimated to be 2 x 10~ 3 M assuming a CO abundance with 
respect to H2 of 10~ 3 (see Tabled- 
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Fig. 9. AARP of the APEX CO(/ = 3 - 2) intensity averaged 
over the central 10km s _1 (squares with errorbars). The model 
CO(7 = 3-2) AARP, where the shell position coincides, within 
the uncertainties, with that of shell 3, is also shown (solid line). 



4.8. The mass of the dust shell 

All light scattered by dust is likely to be optically thin (see be- 
low), and it is therefore possible to do a simple analysis of the 
observational results, based on an analytical approach described 
in GD2001. The ratio of the circumstellar flux over the stellar 
flux is given by 



a eS A/l sc 7V sc cr sc ff e ff A/l sc 



AA f 



4< h 



AA f 



(12) 



where a e ff is ratio of the average stellar flux density in the scatter- 
ing wavelength range to the average stellar flux density over the 
filter width, AA SC and AAf the scattering wavelength range and 
filter width, respectively, N sc the number of scatterers, <x sc the 
scattering cross section, and R s h the radius of the shell. For dust 
scattering the maximum F sc /F+ ratio is equal to one (a e ff = 1, 
and AA SC = AAf), and the measured values are all much lower 
than one (Tabled, hence the scattering is optically thin (assum- 
ing a homogeneous medium). 

Assuming single-sized, spherical dust grains the number of 
scatterers and the scattering cross section are given by 



_3Md_ 

47rp g fl 3 

Que™ 2 , 



(13) 
(14) 



where is the dust mass in the shell, p„ the density of a single 
dust grain, a the radius of a spherical grain, and Q sc the grain 
scattering efficiency. Q sc is determined using Mie scattering the- 
ory. For simplicity isotropic scattering is assumed. Amorphous 
carbon dust grains with optical constants from Suh (2000) were 
adopted. Typical values for dust in the circumstellar envelopes 
of AGB carbon stars are a = 0. 1 fim and p g = 2 g cm -3 . 

The F sc /F+ ratios are given by the observations. The 
model brightness distributions do not reproduce the observed 
I S h AARPs well, likely due to the simplified assumption of a 
constant grain size, the dust density distribution, and a clumpy 



medium. Hence, using the model brightness distributions to de- 
termine F sc /F+ ratios and dust masses individually for shell 3 
and 4 would result in very uncertain results. However, shell 4 
dominates the contribution from dust (Sects. [4~2l and |4~4b . We 
therefore assume that all dust-scattered light comes from shell 
4, and the total observed fluxes are used in determining the dust 
mass (they are corrected for line scattering in the F59 and F77 
filters). The resulting dust mass of shell 4 (averaged over all four 
filters) is 5x 10~ 5 M Q , Tabled Schoier et al. ( 120031 ) modelled the 
thermal emission from U Ant type, and derived a dust mass of 
(1.3 ± 1.2)x 10~ 4 M o ataradiusof (5.1 ± 4.0) x 10 17 cm, con- 
sistent with the mass and radius derived for shell 4 here within 
the considerable uncertainties. Izumiura et al. ([T997) modelled 
two dust shells observed in high-resolution IRAS images and de- 
rived a dust mass for the inner shell of 2.2 x 10~ 5 M Q at a distance 
of 2.0 x 10 17 cm, consistent with the results obtained for shell 4 
here. GD2003 derived a dust mass for shell 4 of 4 x 10~ 6 M Q , a 
factor of 10 lower than our estimate, emphasising the uncertainty 
in estimated dust mass. 



4.9. The mass of the gas shell 

Using the F sc /F* ratios for scattering in the NaD and KI lines, 
corrected for the contribution from dust scattering, it is in princi- 
ple possible to estimate the gas mass of shell 3 based on the an- 
alytical approach from GD2001 (assuming that the contribution 
to the observed total fluxes from the gas is entirely due to shell 
3). The requirement of optically thin scattering sets upper limits 
also to the F^fF* ratios for the KI and NaD lines. For line scat- 
tering we assume the scattering wavelength range to be set by the 
turbulent velocity width. A value typical for circumstellar envi- 
ronments is lkms . For the strong resonance lines a e ff«i0.5 
and the filter widths of the F59 and F77 filters is AAf = 5 nm, 
putting the upper limit to the F sc /F+ ratio for the KI and NaD 
lines at *2.5xl0~ 4 and *4xl0~ 4 , respectively (Equ. [HI The 
observed F sc /F+ ratio is above this limit by a factor of »20 in 
the F59 filter, while the F sc /F+ ratio in the F77 filter lies close 
to the limit (only a factor of « 3 higher). We conclude that scat- 
tering in the NaD lines is optically thick, while scattering in the 
KI lines may be at least partially optically thin. This is a likely 
explanation for the different AARPs in the two filters, as well as 
for the difference in the appearance in the EMMI data. Hence, it 
is reasonable to expect that we can get at most a lower limit to 
the shell 3 gas mass using the optically thin approximation for 
the line-scattered light. 

The best gas mass estimate of shell 3 is consequently ob- 
tained from the CO line modelling, 2 x 10~ 3 M Q . Assuming that 
shells 3 and 4 were created at the same time (see below), the de- 
rived dust and gas masses give a dust-to-gas mass ratio of « 0.02 
for the ejected material, not unreasonable for carbon-rich AGB 
stars. 

We can use the optically thin approximation and the ob- 
served Fsc/F* ratios (corrected for the dust contribution) to 
check whether the Na and K results are consistent with the esti- 
mated shell gas mass. Assuming a Boltzmann distribution of the 
population levels we have for the number of scatterers, 



iV« = 



yum H Z 

and for the line scattering cross section, 

= J 1 q 2 A 2 f 

AA SC 4e m e c 2 ' 



(15) 



(16) 
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Table 5. Derived shell masses based on the EFOSC2 data (dust) 
and the APEX CO line data (gas). 
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where r/ n is the fraction of neutral atoms, fx the fractional abun- 
dance of species X relative to H, g\ the statistical weight of the 
lower level, E\ the energy of the lower level, Z the partition func- 
tion, q the charge of an electron, / the oscillator strength of the 
line, and the other parameters follow usual notations. Assuming 
solar abundances of Na and K we find that only 13% and 27% of 
the Na and K, respectively, is required to remain in neutral form 
to explain the observational results. 

5. Discussion and conclusions 

5.1. The identification of shells 

Based on the EFOSC2 data we clearly detect two geometrically 
thin shells around U Ant with radii of 43" and 50" (shells 3 and 
4, respectively). Although shell 4 is relatively weak in the im- 
ages of total scattered light, it is clearly visible in the images of 
polarised light. We also confirm the tentative detection of addi- 
tional structure inside shell 3 in the F59 filter images, a result 
supported by the EMMI NaD data. We proceed to discuss the 
tentative shells 1 and 2, and the shells 3 and 4 in separate sec- 
tions below. 

5.1.1. Shells 1 and 2 

These shells were tentatively introduced by GD2003. The struc- 
tures attributed to these shells were only seen in the EFOSC2 
F59 data. The I s y, image shows signs of arcs that may form an 
essentially circular shell (shell 2). However, this putative shell is 
not centered on the star. There are also arcs that seem to 'con- 
nect' this shell and shell 3. The presence of shell 1 was inferred 
only from the /,/, AARP. However, here we provide evidence for 
structures in the radial ranges of shells 1 and 2 also in the EMMI 
NaD line data, but these structures do not form complete (or even 
partial) ellipses as expected from spherical shells centered on the 
star. Curiously, these structures are not seen in any of the other 
data. The most reasonable explanation for this is that the optical 
depth of the NaD line scattering is high (above one) and that of 
the other scattering agents is low (below one). This is consistent 
with the asymmetry in the brightness between the red- and blue- 
shifted sides in the EMMI data. The effect is stronger in the NaD 
lines, indicating a higher optical depth than in the KI line. The 
spatial resolution and the dynamic range of the CO data are not 
enough to identify structures in this radial range. We conclude 
that there is no strong evidence for the existence of full shells 
inside shell 3, but structures do exist and their interpretation is 
not clear. Whether these structures indeed are closer to the star, 
or whether they are structures within the observed shells is also 
not clear. 

5.1.2. Shells 3 and 4 

Shells 3 and 4 are clearly present as separate entities in the po- 
larisation data at radii of » 43" and 50", respectively. In the P s h 
images shells 3 and 4 have a large degree of circular symme- 



try, with slight variations in the intensity indicating a somewhat 
clumpy structure. There is a gap in shells 3 and 4 in the north- 
eastern corner of the polarised images, however, this is likely 
due to artefacts in the images, rather than a gap in the density 
distribution in the shells. The single-dish CO radio observations 
do not resolve individual shells, nor do they indicate a multiple- 
shell structure. However, there are several indications that the 
CO line emission comes predominantly from shell 3 (Sect. l4~7l >. 
In addition, there is a 'hole' in the south-western quadrant in the 
CO maps (Olofsson et al. |1996b and this is seen also in the F59 
and F77 images. Based on the flux densities in the different fil- 
ters and the polarisation data, it is clear that shell 4 essentially 
consists of only dust, while shell 3 is dominated by gas. 

The results of the brightness distribution fits give FWHM 
widths (of the density distributions of « 2" and » 7" for shell 3 
and 4, respectively. The narrow width of shell 3 is confirmed by 
the EMMI data, and consistent with the CO line data. Thus, we 
conclude that the gas shell is geometrically thin (AR/R<0.1), 
while the external dust shell is wider. 

5.2. The evolution of the detached shells 

Assuming a distance of 260 pc to U Ant (the Hipparcos dis- 
tance), we derive mean shell sizes of 1.7 x 10 17 cm and 1.9 x 
10 17 cm for shells 3 and 4, respectively. The EMMI data shows 
that shell 3 has an expansion velocity of « 19.5 km s -1 , a value 
consistent with the CO line results. This corresponds to a dy- 
namical shell 3 age of as 2700 years, assuming that the shell ex- 
pansion velocity has been constant with time. 

We believe shell 4 to consist of only dust, while shell 3 is 
dominated by gas. A possible interpretation of this fact is that 
both shells originate in the same event, but that the dust, mov- 
ing at a slightly higher velocity, has to a large extent separated 
from the gas and formed a shell of its own further out. This was 
already suggested by GD2003. Assuming that both shells were 
created 2700 years ago, this implies a drift velocity of 2.8 km s _1 
between the gas and the dust, certainly reasonable for winds 
of AGB stars with high mass-loss rates (Ramstedt et al. 12008b . 
Further, based on the results for the gas mass of shell 3 from the 
CO line models, 2 x 10~ 3 M Q , and assuming that all the dust in 
shells 3 and 4 comes from the same event, the dust-to-gas mass 
ratio of the ejected material is « 0.02. 

The derivation of the drift velocity assumes that the dust, 
once separated from the gas shell, moves at a constant velocity. 
However, when the dust leaves the gas, the drag force due to col- 
lisions with the gas disappears, and the radiation pressure from 
the star again accelerates the grains. The equation of motion for 
grains that have separated from the gas is 

tL = ^LQsc 
dt 2 16nap g cr 2 

where r is the distance from the star, L the stellar luminosity, Q sc , 
a, and p g the scattering efficiency, the radius, and the density of a 
grain, respectively. We solve Equ.[17]numerically, assuming the 
same grain properties as in Sect. 14.81 to determine at what time 
shell 4 must have separated from shell 3 in order for both shells 
to arrive at their observed distances from the star. For the evo- 
lution of the shells after separation, we assume that shell 3 con- 
tinues to expand with a constant velocity of « 19.5 km s _1 , while 
shell 4 has a velocity of 22.3kms _I at the moment of separa- 
tion. The solution shows that shell 4 and 3 must have separated 
only 1 10 years ago for both shells to reach their current posi- 
tions. However, the probability of observing the detached dust 
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and gas shells immediately after separation is unlikely, and it is 
probable that the dust separated at an earlier time. Therefore, we 
would expect the dust shell to be at a significantly larger distance 
than what is observed. There could be several reasons why the 
distance between shell 3 and 4 is smaller than expected. A possi- 
ble explanation is friction outside the gas shell, due to an earlier 
stellar wind. Simple estimates suggest, however, that this is not 
efficient enough. Another possibility is that the grain radius a is 
greater than the adopted value of 0. 1 fim. Since the acceleration 
scales as a -1 , the predicted distance between the shells would 
shrink. The large width of the dust shell is possibly attributable 
to dust grains of different sizes moving with different velocities 
as well. Finally, our assumption concerning the shells, having 
constant mass, constant thickness, and a smooth medium, clearly 
is an oversimplification, as shown by the hydrodynamic simula- 
tions by Steffen et al. (fT998] >, Steffen & Schonberner (120001 ), and 
Mattsson et al. (120071 ). The gradual sweeping up of earlier stel- 
lar winds may increase the mass density in the gas shell and thus 
the friction with time (Schoier et al. |20051 l, which could lead to a 
slower passage of the dust grains through the shell, hence delay- 
ing the separation of dust and gas. A realistic modelling of the 
dust drift in such a system is beyond the scope of this paper. 



5.3. The origin of the detached shells 

A link between geometrically thin, detached shells and thermal 
pulses was proposed already 20 years ago, although the exact 
mechanism remained unclear (Olofsson et al. 1 19901 1. The shells 
were suggested to have formed during a period of increased 
mass-loss rate due to changes in the stellar temperature and lu- 
minosity during the thermal pulse. However, hydrodynamical 
models showed that a mass injection due to a period of high 
mass loss alone proved not sufficient to form the observed gas 
shells (Steffen et al. 119981 ), as they tended to diffuse. Instead, 
the interaction of the thermal pulse wind with a previous wind 
of lower velocity proved to be a more promising mechanism for 
producing the shells (Steffen & Schonberger |2000l l. Although 
Steffen & Schonberger managed to create detached shells in 
connection with thermal pulses in their models, they relied on 
empirical mass-loss descriptions, even during the thermal pulse. 
Mattson et al. (2007 ) for the first time self-consistently modelled 
the change in stellar parameters, the dynamics in the wind accel- 
eration zone, and the evolution of the expanding circumstellar 
envelope during a thermal pulse. Their models form detached 
shells due to an increase in expansion velocity and mass loss 
during the pulse. Critical for the formation of detached shells in 
their models is the effect of the stellar pulsation characteristics 
(described as a piston at as 0.9 stellar radii), and the amplitude 
of the change in mass-loss rate and expansion velocity. They 
find shell densities of the order of 10~ 21 gem 3 at distances of 
1 -2xl0 17 cm, and expansion velocities of 18kms~', in good 
agreement with the values derived from our observations. The 
resulting relative shell thickness in the model (AR/R « 0.01) 
is approximately a factor of 10 smaller than what observations 
show. However, the thickness of a shell is sensitive to model as- 
sumptions, such as heating and cooling, isothermal or adiabatic 
shocks, and the jump in expansion velocity. Observational sup- 
port for a two-wind interaction scenario was given by radiative 
transfer modelling of dust and CO radio line emission for seven 
stars with geometrically thin shells by Schoier et al. ( 120051 ). They 
found evidence of increasing shell mass and decreasing shell ex- 
pansion velocities with increasing shell radius. This is consis- 
tent with a scenario where an expanding shell collides with a 



previous, slower wind, sweeping up material and hence slowing 
down. 

Although a separation of the dust and gas due to different 
expansion velocities seems a reasonable explanation, a different 
behaviour of the dust and gas components may also be explained 
by different creation mechanisms. Steffen et al. (119981 ) describe 
the creation of dust shells as a natural consequence of the large 
drop in mass-loss rate after the thermal pulse. The decrease in 
mass-loss rate leads to a sharp decrease of dust condensation 
during the mass-loss rate minimum after the thermal pulse, in 
turn leading to an uncoupling of the dust and gas. This increases 
the dust outflow velocity resulting in the formation of a detached 
dust shell, i.e. interaction with a previous wind is not necessary 
to create detached shells of dust. 

We conclude that the results from the observations are con- 
sistent with a thermal pulse scenario (as described above) taking 
place about 2700 years ago. A gas/dust shell was created that 
expanded with a velocity higher than during normal mass-loss 
conditions. The dust moved with a slightly higher velocity and 
it remained unaffected by any dust in a previous wind. The gas, 
on the other hand, interacted with the previous lower velocity 
wind creating a narrow shell with gradually increasing gas mass. 
The question arises whether the gas or the dust density distribu- 
tion follows most closely the mass-loss rate modulation during 
the ejection. The gas shell width corresponds to a time scale of 
about 150 years, while the dust shell width corresponds to about 
400 years. The gas shell is affected by interaction while the dust 
shell can be affected by grain-size-dependent expansion veloci- 
ties. To answer this question requires models that follow in detail 
the evolution of the gaseous and dusty media. 

We confirm that there is circumstellar structure in the re- 
gion inside shell 3 (i.e., at projected distances closer to the star), 
but there is no good evidence that these form shells, i.e., that 
they can be attributed to additional variations in the mass-loss 
rate. Higher quality data is required to solve this issue, which 
remains a point of some concern for the overall interpretation. 
Assuming the same expansion velocity for the putative shell 2 as 
for shell 3 implies an age difference of only 600 years - not in 
any obvious way consolable with the time scale between consec- 
utive thermal pulses. Multiple shells (or rather arcs) have been 
observed around e.g. the carbon star IRC+10216 (Mauron & 
Huggins [1999l Mauron & Hug gins |2000l Fong et al. |2003l , sep- 
arated by time scales of a few hundred to thousand years. The 
creation of this type of shells with time scales of less than the 
interpulse time scale may be explained by density and velocity 
modulations during a single thermal pulse, in addition to being 
the consequence of successive pulses (Villaver et al. 120021 ). 
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